Azimuthal anisotropy of n° production in Au+Au collisions at */5nn = 200 GeV: 
Path-length dependence of jet-quenching and the role of initial geometry 



A. Adare, 11 S. Afanasiev, 26 C. Aidala, 39 N.N. Ajitanand, 56 Y. Akiba, 50 - 51 H. Al-Bataineh, 45 J. Alexander, 56 
K. Aoki, 32 * 50 Y. Aramaki, 10 E.T. Atomssa, 33 R. Averbeck, 57 T.C. Awes, 46 B. Azmoun, 5 V. Babintsev, 22 
M. Bai, 4 G. Baksay, 18 L. Baksay, 18 K.N. Barish, 6 B. Bassalleck, 44 A.T. Basye, 1 S. Bathe, 6 V. Baublis, 49 
C. Baumann, 40 A. Bazilevsky, 5 S. Belikov, 5 'B R. Belmont, 61 R. Bennett, 57 A. Berdnikov, 53 Y. Bcrdnikov, 53 
A. A. Bickley, 11 J.S. Bok, 64 K. Boyle, 57 M.L. Brooks, 35 H. Buesching, 5 V. Bumazhnov, 22 G. Bunce, 5 ' 51 S. Butsyk, 35 
CM. Camacho, 35 S. Campbell, 57 C.-H. Chen, 57 CY. Chi, 12 M. Chiu, 5 I.J. Choi, 64 R.K. Choudhury, 3 
P. Christiansen, 37 T. Chujo, 60 P. Chung, 56 O. Chvala, 6 V. Cianciolo, 46 Z. Citron, 57 B.A. Cole, 12 M. Connors, 57 
P. Constantin, 35 M. Csanad, 16 T. Csorgo, 29 T. Dahms, 57 S. Dairaku, 32 ' 50 I. Danchev, 61 K. Das, 19 A. Datta, 39 
G. David, 5 A. Denisov, 22 A. Deshpande, 51 ' 57 E.J. Desmond, 5 O. Dietzsch, 54 A. Dion, 57 M. Donadelli, 54 
O. Drapier, 33 A. Drees, 57 K.A. Drees, 4 J.M. Durham, 57 A. Durum, 22 D. Dutta, 3 S. Edwards, 19 Y.V. Efremenko, 46 

F. Ellinghaus, 11 T. Engelmore, 12 A. Enokizono, 34 H. En'yo, 50 ' 51 S. Esumi, 60 B. Fadem, 41 D.E. Fields, 44 
M. Finger, Jr., 7 M. Finger, 7 F. Fleuret, 33 S.L. Fokin, 31 Z. Fraenkel, 63 'B J-E. Frantz, 57 A. Franz, 5 A.D. Frawley, 19 
K. Fujiwara, 50 Y. Fukao, 50 T. Fusayasu, 43 L Garishvili, 58 A. Glenn, 11 H. Gong, 57 M. Gonin, 33 Y. Goto, 50 ' 51 
R. Granier de Cassagnac, 33 N. Grau, 12 S.V. Greene, 61 M. Grosse Perdekamp, 23 - 51 T. Gunji, 10 H.-A. Gustafsson, 37 'B 
J.S. Haggerty, 5 K.I. Hahn, 17 H. Hamagaki, 10 J. Hamblen, 58 J. Hanks, 12 R. Han 48 E.P. Hartouni, 34 E. Haslum, 37 
R. Hayano, 10 M. Heffner, 34 S. Hegyi, 29 T.K. Hemmick, 57 T. Hester, 6 X. He, 20 J.C. Hill, 25 M. Hohlmann, 18 
W. Holzmann, 12 K. Homma, 21 B. Hong, 30 T. Horaguchi, 21 D. Hornback, 58 S. Huang, 61 T. Ichihara, 50 ' 51 
R. Ichimiya, 50 J. Ide, 41 Y. Ikeda, 60 K. Imai, 32 ' 50 M. Inaba, 60 D. Isenhower, 1 M. Ishihara, 50 T. Isobe, 10 M. Issah, 61 
A. Isupov, 26 D. Ivanischev, 49 B.V. Jacak, 57 -Q J. Jia, 5 ' 56 J. Jin, 12 B.M. Johnson, 5 K.S. Joo, 42 D. Jouan, 47 
D.S. Jumper, 1 F. Kajihara, 10 S. Kametani, 50 N. Kamihara, 51 J. Kamin, 57 J.H. Kang, 64 J. Kapustinsky, 35 

D. Kawall, 39 ' 51 M. Kawashima, 52 ' 50 A.V. Kazantsev, 31 T. Kempel, 25 A. Khanzadeev, 49 K.M. Kijima, 21 B.I. Kim, 30 
D.H. Kim, 42 D.J. Kim, 27 E.J. Kim, 8 E. Kim, 55 S.H. Kim, 64 Y.J. Kim, 23 E. Kinney, 11 K. Kiriluk, 11 A. Kiss, 16 

E. Kistenev, 5 L. Kochenda, 49 B. Komkov 49 M. Konno, 60 J. Koster, 23 D. Kotchetkov 44 A. Kozlov, 63 A. Krai, 13 
A. Kravitz, 12 G.J. Kunde, 35 K. Kurita, 52 * 50 M. Kurosawa, 50 Y. Kwon, 64 G.S. Kyle, 45 R. Lacey, 56 Y.S. Lai, 12 
J.G. Lajoie, 25 A. Lebedev, 25 D.M. Lee, 35 J. Lee, 17 K.B. Lee, 30 K. Lee, 55 K.S. Lee, 30 M.J. Leitch, 35 M.A.L. Leite, 54 
E. Leitner, 61 B. Lenzi, 54 P. Liebing, 51 L.A. Linden Levy, 11 T. Liska, 13 A. Litvinenko, 26 H. Liu, 35 ' 45 M.X. Liu, 35 
X. Li, 9 B. Love, 61 R. Luechtenborg 40 D. Lynch, 5 C.F. Maguire, 61 Y.I. Makdisi 4 A. Malakhov, 26 M.D. Malik, 44 
V.I. Manko, 31 E. Mannel, 12 Y. Mao, 48 - 50 H. Masui, 60 F. Matathias, 12 M. McCumber, 57 PL. McGaughey, 35 
N. Means, 57 B. Meredith, 23 Y. Miake, 60 A.C. Mignerey, 38 P. Mikes, 7 ' 24 K. Miki, 60 A. Milov, 5 M. Mishra, 2 
J.T. Mitchell, 5 A.K. Mohanty, 3 Y. Morino, 10 A. Morreale, 6 D.P Morrison, 5 T.V. Moukhanova, 31 J. Murata, 52 ' 50 

S. Nagamiya, 28 J.L. Nagle, 11 M. Naglis, 63 M.I. Nagy, 16 I. Nakagawa, 50 ' 51 Y. Nakamiya, 21 T. Nakamura, 21 * 28 
K. Nakano, 50 ' 59 J. Newby, 34 M. Nguyen, 57 R. Nouicer, 5 A.S. Nyanin, 31 E. O'Brien, 5 S.X. Oda, 10 CA. Ogilvie, 25 
K. Okada, 51 M. Oka, 60 Y. Onuki, 50 A. Oskarsson, 37 M. Ouchida, 21 K. Ozawa, 10 R. Pak, 5 V. Pantuev, 57 
V. Papavassiliou, 45 I.H. Park, 17 J. Park, 55 S.K. Park, 30 W.J. Park, 30 S.F. Pate, 45 H. Pei, 25 J.-C. Peng, 23 
H. Pereira, 14 V. Peresedov, 26 DYu. Peressounko, 31 C. Pinkenburg, 5 R.P. Pisani, 5 M. Proissl, 57 M.L. Purschke, 5 

A.K. Purwar, 35 H. Qu, 20 J. Rak, 27 A. Rakotozafindrabe, 33 I. Ravinovich, 63 K.F. Read, 46 ' 58 K. Reygers, 40 
V. Riabov, 49 Y. Riabov, 49 E. Richardson, 38 D. Roach, 61 G. Roche, 36 S.D. Rolnick, 6 M. Rosati, 25 CA. Rosen, 11 
S.S.E. Rosendahl, 37 P. Rosnet, 36 P. Rukoyatkin, 26 P. Ruzicka, 24 B. Sahlmueller, 40 N. Saito, 28 T. Sakaguchi, 5 
K. Sakashita, 50 ' 59 V. Samsonov, 49 S. Sano, 10 * 62 T. Sato, 60 S. Sawada, 28 K. Sedgwick, 6 J. Seele, 11 R. Seidl, 23 

AYu. Semenov, 25 R. Seto, 6 D. Sharma, 63 I. Shein, 22 T.-A. Shibata, 50 ' 59 K. Shigaki, 21 M. Shimomura, 60 
K. Shoji, 32 ' 50 P. Shukla, 3 A. Sickles, 5 C.L. Silva, 54 D. Silvermyr, 46 C. Silvestre, 14 K.S. Sim, 30 B.K. Singh, 2 
CP. Singh, 2 V. Singh, 2 M. Slunecka, 7 R.A. Soltz, 34 W.E. Sondheim, 35 S.P. Sorensen, 58 I.V. Sourikova, 5 
N.A. Sparks, 1 P.W. Stankus, 46 E. Stenlund, 37 S.P. Stoll, 5 T. Sugitate, 21 A. Sukhanov, 5 J. Sziklai, 29 

E. M. Takagui, 54 A. Taketani, 50 ' 51 R. Tanabe, 60 Y. Tanaka, 43 K. Tanida, 32 ' 50 * 51 M.J. Tanncnbaum, 5 S. Tarafdar, 2 

A. Taranenko, 56 P. Tarjan, 15 H. Themann, 57 T.L. Thomas, 44 M. Togawa, 32 ' 50 A. Toia, 57 L. Tomasek, 24 
H. Torii, 21 R.S. Towell, 1 I. Tserruya, 63 Y. Tsuchimoto, 21 C. Vale, 5 ' 25 H. Valle, 61 H.W. van Hecke, 35 
E. Vazquez-Zambrano, 12 A. Veicht, 23 J. Velkovska, 61 R. Vertesi, 15,29 A. A. Vinogradov, 31 M. Virius, 13 
V. Vrba, 24 E. Vznuzdaev, 49 X.R. Wang, 45 D. Watanabe, 21 K. Watanabe, 60 Y. Watanabe, 50 ' 51 F. Wei, 25 



2 



R. Wei, 56 J. Wessels, 40 S.N. White, 5 D. Winter, 12 J.P. Wood, 1 C.L. Woody, 5 R.M. Wright, 1 M. Wysocki, 11 

W. Xic, 51 Y.L. Yamaguchi, 10 K. Yamaura, 21 R. Yang, 23 A. Yanovich, 22 J. Ying, 20 S. Yokkaichi, 50 ' 51 
G.R. Young, 46 I. Younus, 44 Z. You, 48 I.E. Yushmanov, 31 W.A. Zajc, 12 C. Zhang, 46 S. Zhou, 9 and L. Zolin 26 

(PHENIX Collaboration) 

Abilene Christian University, Abilene, Texas 79699, USA 
2 Department of Physics, Banaras Hindu University, Varanasi 221005, India 
Bhabha Atomic Research Centre, Bombay 400 085, India 
4 Collider- Accelerator Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
5 Physics Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
6 University of California - Riverside, Riverside, California 92521, USA 
'Charles University, Ovocny trh 5, Praha 1, 116 36, Prague, Czech Republic 
Chonbuk National University, Jeonju 561-756, Korea 
'China Institute of Atomic Energy (CIAE), Beijing, People's Republic of China 
10 Center for Nuclear Study, Graduate School of Science, University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 

11 University of Colorado, Boulder, Colorado 80309, USA 
12 Columbia University, New York, New York 10027 and Nevis Laboratories, Irvington, New York 10533, USA 
1:1 Czech Technical University, Zikova 4, 166 36 Prague 6, Czech Republic 
14 Dapnia, CEA Saclay, F-91191, Gif-sur-Yvette, France 
1,1 Debrecen University, H-4010 Debrecen, Egyetem ter 1, Hungary 
16 ELTE, Eotvos Lordnd University, H - 1117 Budapest, Pdzmdny P. s. 1/A, Hungary 
17 Ewha Womans University, Seoul 120-750, Korea 
18 Florida Institute of Technology, Melbourne, Florida 32901, USA 
19 Florida State University, Tallahassee, Florida 32306, USA 
20 Georgia State University, Atlanta, Georgia 30303, USA 
21 Hiroshima University, Kagamiyama, Higashi- Hiroshima 739-8526, Japan 
22 IHEP Protvino, State Research Center of Russian Federation, Institute for High Energy Physics, Protvino, 142281, Russia 
23 University of Illinois at Urbana- Champaign, Urbana, Illinois 61801, USA 
24 Institute of Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, 182 21 Prague 8, Czech Republic 

25 Iowa State University, Ames, Iowa 50011, USA 
26 Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia 
27 Helsinki Institute of Physics and University of Jyvdskyld, P.O.Box 35, FI-40014 Jyvdskyld, Finland 
28 KEK, High Energy Accelerator Research Organization, Tsukuba, Ibaraki 305-0801, Japan 
29 KFKI Research Institute for Particle and Nuclear Physics of the Hungarian Academy 
of Sciences (MTA KFKI RMKI), H-1525 Budapest 114, POBox 49, Budapest, Hungary 
"'Korea University, Seoul 136-701, Korea 
31 Russian Research Center "Kurchatov Institute", Moscow, Russia 
32 Kyoto University, Kyoto 606-8502, Japan 
33 Laboratoire Leprince-Ringuet, Ecole Poly technique, CNRS-IN2P3, Route de Saclay, F-91128, Palaiseau, France 
34 Lawrence Livermore National Laboratory, Livermore, California 94550, USA 
35 Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
36 LPC, Universite Blaise Pascal, CNRS-IN2P3, Clermont-Fd, 63177 Aubiere Cedex, France 
37 Department of Physics, Lund University, Box 118, SE-221 00 Lund, Sweden 
38 University of Maryland, College Park, Maryland 20742, USA 
19 Department of Physics, University of Massachusetts, Amherst, Massachusetts 01003-9337, USA 
41 ' Institut fur Kernphysik, University of Muenster, D-48149 Muenster, Germany 
41 Muhlenberg College, Allentown, Pennsylvania 18104-5586, USA 
42 Myongji University, Yongin, Kyonggido 449-728, Korea 
43 Nagasaki Institute of Applied Science, Nagasaki-shi, Nagasaki 851-0193, Japan 
44 University of New Mexico, Albuquerque, New Mexico 87131, USA 
45 New Mexico State University, Las Cruces, New Mexico 88003, USA 
4< 'Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 
47 IPN-Orsay, Universite Paris Sud, CNRS-IN2P3, BP1, F-91406, Orsay, France 
48 Peking University, Beijing, People's Republic of China 
49 PNPI, Petersburg Nuclear Physics Institute, Gatchina, Leningrad region, 188300, Russia 
5 "RIKEN Nishma Center for Accelerator-Based Science, Wako, Saitama 351-0198, JAPAN 
51 RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
52 Physics Department, Rikkyo University, 3-34-1 Nishi-Ikebukuro, Toshima, Tokyo 171-8501, Japan 
53 Saint Petersburg State Polytechnic University, St. Petersburg, Russia 
54 Universidade de Sao Paulo, Instituto de Fisica, Caixa Postal 66318, Sao Paulo CEP05315-970, Brazil 
55 Seoul National University, Seoul 151-742, Korea 
S6 Chemistry Department, Stony Brook University, Stony Brook, SUNY, New York 11794-3400, USA 
57 Department of Physics and Astronomy, Stony Brook University, SUNY, Stony Brook, New York 11794, USA 



3 



University of Tennessee, Knoxville, Tennessee 37996, USA 
' 9 Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro, Tokyo 152-8551, Japan 
60 Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
61 Vanderbilt University, Nashville, Tennessee 37235, USA 
62 Waseda University, Advanced Research Institute for Science and 
Engineering, 17 Kikui-cho, Shinjuku-ku, Tokyo 162-0044, Japan 
63 Weizmann Institute, Rehovot 76100, Israel 
64 Yonsei University, IPAP, Seoul 120-749, Korea 
(Dated: June 21, 2010) 

We have measured the azimuthal anisotropy of 7v° production for 1 < pt < 18 GeV/c for Au+Au 
collisions at y / %N = 200 GeV. The observed anisotropy shows a gradual decrease for 3 < pr < 7— 10 
GeV/c, but remains positive beyond 10 GeV/c. The magnitude of this anisotropy is under-predicted, 
up to at least ~ 10 GeV/c, by current perturbative QCD (pQCD) energy-loss model calculations. 
An estimate of the increase in anisotropy expected from initial-geometry modification due to gluon 
saturation effects and fluctuations is insufficient to account for this discrepancy. Calculations that 
implement a path length dependence steeper than what is implied by current pQCD energy-loss 
models show reasonable agreement with the data. 

PACS numbers: 25.75.Dw 



A central goal of high-energy nuclear physics is to un- 
derstand the properties of the strongly-coupled Quark 
Gluon Plasma (sQGP), a new form of nuclear mat- 
ter identified at the Relativistic Heavy Ion Collider 
(RHIC) [1] . An important tool for this goal is jet quench- 
ing or the suppression of high transverse momentum (px) 
hadron yields as a result of in-medium energy loss of high 
Px jets [2J. Such suppression was first observed in mea- 
surements of the nuclear modification factor for single 
hadron yields Raa = /T rfjv A A — , where dN^A is the dif- 



ferential yield in Au+Au collisions, da pp is the differential 
cross section in p+p collisions for given px, and (Taa) is 
the nuclear overlap integral for given Au+Au centrality 
bin Q. Later on this effect was also observed in mea- 
surements of dihadron Q and 7-hadron correlations [f| . 

Current theoretical descriptions of jet quenching are 
commonly based on a pertubative QCD (pQCD) frame- 
work 6), which assumes that the coupling of jets with 
the medium is weak, even though the medium itself 
is strongly-coupled. Prompted by the large amount 
of experimental data from RHIC, several sophisticated 
pQCD-based models have been developed in the last 
decade @, [f|. These models have provided initial esti- 
mates of the properties of the sQGP, such as the mo- 
mentum broadening per mean free path, d = (k T )/X, 
and the energy loss per unit length, dE/dl [a-lg]. 

Despite these successes, the pQCD description of jet 
quenching faces several challenges (see Ref @). Be- 
sides a large discrepancy among models of extracted 
medium properties such as q [3], the energy-loss mod- 
els also disagree in their predictions of the azimuthal 
anisotropy of high px hadrons [H . The latter character- 
izes hadron emission relative to the reaction plane angle 
(*hp), dN/d(<P-V R p) cx (l + 2w 2 cos(2(0-* R p))). Such 
azimuthal anisotropy ensues because the hadron yield 
is more suppressed along the long axis of the almond- 
shaped fireball than the short axis. Thus the magnitude 



of the anisotropy, i>2, is sensitive to the path length (/) 
dependence of energy loss, which scales as AE ~ I for col- 
lisional energy loss [lCf , AE ~ I 2 for coherent radiative 
energy loss [lCf , and AE ~ I 3 for a non-perturbative en- 
ergy loss calculation using AdS/CFT gravity-gauge dual 
theory However, our ability to probe such / depen- 
dences hinges not only on precision data at high px, but 
also on a good understanding of the role of the initial col- 
lision geometry. One geometry commonly used in energy- 
loss models is based on the Optical Glauber model [12j . 
which assumes a smooth Woods-Saxon nuclear geometry. 
Such geometry ignores the event-by-event shape distor- 
tion due to spacial fluctuations of participating nucle- 
ons [l3j], and a possible overall shape distortion due to 
gluon saturation effects, i.e. the so called CGC geome- 
try [13]. These effects have been shown to be important 
(up to 15-30% each) for elliptic flow at low px [HI fl6j - 
However, their influences on jet quenching V2 at high px 
are not well studied to date. 

In this Letter we present a new measurement of the ir° 
anisotropy in ^/snn = 200 GeV Au+Au collisions. This 
measurement complements our prior results [IMl, but 
significantly increases both the px reach and the statis- 
tical precision above 6 GeV/c, allowing for quantitative 
comparisons to energy-loss models, as well as detailed 
investigations of the role of the initial collision geometry. 

Results were obtained from ~ 3.5 x 10 9 minimum bias 
events taken in 2007. Event centrality was determined 
by the number of charged particles detected in the Beam- 
Beam Counters (BBGJB.O < \r]\ < 3.9). A Monte-Carlo 
(MC) Glauber model [12j was used to estimate the aver- 
age number of participating nucleons (A par t) and (Iaa) 
for each centrality class. 

Previous PHENIX analyses [l9[ estimated the RP us- 
ing the charged particles detected in the BBC. Several 
new detectors, installed symmetrically on both sides of 
the beam line, provided additional RP measurements in 
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FIG. 1: (a)-(f): tt° V2 using reaction plane determined with MPC and RXNi n combined as a function of pt for different 
centralities. (g)-(i): ratios of V2 measured separately using MPC (solid triangles) and RXNi n (open triangles) to combined 
result; the dashed lines indicate the systematic error. 



2007: the Muon Piston Calorimeters (MPC, 3.1 < \r}\ < 
3.9) and the Reaction Plane detectors in two r\ ranges, 
RXN in (RXN out ) in 1.5(1.0) < \r\\ < 2.8(1.5). Each MPC 
is equipped with PDWO4 crystal scintillators to detect 
both charged and neutral particles. Each RXN consists 
of 12 azimuthally segmented paddle scintillators. This 
analysis estimates the RP angle using both the MPC 
and RXN in to provide good resolution, while minimiz- 
ing the potential biases from jets and dijets [20} . The 
error on the RP angle A>i?, and the RP dispersion fac- 
tor orp = (cos 2 A*) are estimated by the sub-event 
method [l9j], giving erp,p ~ 0.52 and 0.73 in the central 
and mid-central collisions, respectively, which is ~ 80% 
better than that for the BBCs. The large dataset and 
improved orp give an equivalent of ~ 15 fold increase in 
statistics over the previous measurement of V2 fl9j ] . 

The methodology for V2 extraction follows our previ- 
ous work (ljj . We reconstruct the neutral pions via the 
7T° — > 7 + 7 decay channel with photons detected in the 
Electromagnetic Calorimeter (EMC, \rj\ < 0.35). We ap- 
ply shower shape and pair asymmetry cuts to reduce the 
combinatorial background. The remaining background 
is subtracted by the mixed event method [19J. The az- 
imuthal distribution of the tt° yields relative to the esti- 
mated RP angle, A(j) — <j) — ^rp, is divided into 6 bins in 
the interval of [0, tt/2], and fit to A^ (l + 2^ aw cos(2A0)). 
The V2 is then obtained by applying the dispersion cor- 
rection V2 = V2 aw /a"BP for each centrality and pr selec- 
tion. The main sources of systematic uncertainties come 
from ctrp and v r j™ . The former is estimated by com- 
paring measurements from different RP detectors, giving 
~ 10% for central and peripheral collisions and ~ 5% 
for mid-central collisions. The latter accounts for the de- 
pendence of V2 on 7T° identification cuts, as well as the 
variation among different sectors of EMC and different 
run groups, and is correlated in px', it is estimated to be 
10% for central collisions and 3% for other collisions. 



Figured] (a)-(f) shows v 2 (pt) for six centrality bins, 
spanning 1-18 GeV/c. In the 10-50% centrality range, 
where the signal is large and the uncertainty is small, 
the i»2 values above 3 GeV/c indicate a slow decrease up 
to 7-10 GeV/c, and remain significantly above zero at 
higher px- The ratios in Fig.[T](g)-(i) confirm the consis- 
tency of i>2 measured using the RP from the MPC or the 
RXN; n and imply that the influence of rapidity depen- 
dent jet bias to the RP, if any, is within the statistical or 
systematic uncertainty of the measurement. 

Figure [2] (a)- (b) shows the centrality dependence of V2 
in two high pt selections. They are compared with four 
pQCD jet quenching model calculations, AMY, HT and 
ASW from [1] and WHDG from [2% The WHDG model 
was calculated for gluon density dN/dg = 1000 — 1600, 
a range constrained by 0-5% (N part = 351) tt° Raa 
data [7j. The calculation assumes analytical Woods- 
Saxon nuclear geometry with a longitudinal Bjorken ex- 
pansion. The AMY, HT and ASW models were fitted 
independently to the 0-5% ir° Raa data [U; they were 
implemented in a 3D ideal hydrodynamic code with iden- 
tical initial Wood-Saxon nuclear geometry, medium evo- 
lution and fragmentation functions. The HT and ASW 
models include only coherent radiative energy loss, while 
the AMY and WHDG also include collisional energy loss. 
The ASW and WHDG models predict quite sizable, and 
similar V2, while the HT and AMY models tend to give 
much smaller V2- However, all models significantly under- 
predict the V2 data in 6 < pr < 9 GeV/c range. For 
Pt > 9 GeV/c, ASW and WHDG results show a bet- 
ter agreement with the 20-30% (A part = 167) centrality 
bin due to a slow decrease of «2 with px in this bin (see 
Fig. QJb)). However, this seems to be accidental, since 
the i>2 values for the other centrality bins remain large, 
and are significantly above the WHDG calculations (the 
p- value for the agreement is < 10 -4 ). 

In all these models, the inclusive suppression Raa and 
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FIG. 2: (Color online) (a)-(b) V2 vs Afpart in two pr ranges; 
(c)-(d) Raa vs ^part in same pr ranges. Each are compared 
with four pQCD models from [| (AMY, HT, ASW) and [H 
(WHDG). Log-scale is used for Raa to better visualize var- 
ious model calculations. Note that the 4r- = 1000 of WHDG 

dg 

corresponds to lower (upper) boundary of the shaded bands 
for V2 (Raa), while 7^ = 1600 corresponds to upper (lower) 
boundary for V2 (Raa)- 



V2 are anti-correlated, i.e. a, smaller Raa implies a larger 
V2 and vice versa. Consequently, more information can 
be obtained by comparing the data with a given model for 
both Raa and v%. Fig. |2](c)-(d) compares the centrality 
dependence of ir° Raa data to four model calculations 
for the same two pr ranges [2l[. The calculations are 
available for a broad centrality range for WHDG, but 
only in 0-5% and 20-30% centrality bins for AMY, HT 
and ASW. The level of agreement varies among the mod- 
els. The HT calculations are slightly above the data in 
the most central bin, while WHDG systematically under- 
predicts the data over the full centrality range, though 
better agrement with the data is obtained for p? > 9 
GeV/c. On the other hand, ASW and AMY calcula- 
tions agree with the data very well in both pt ranges. 
The different levels of agreement among the models are 
partially due to their different trends of Raa with px'. 
WHDG and ASW results have stronger pr dependences 
than what is impled by the data, and tend to deviate 
at low pt when fitted to the full pr range @, HI ■ Given 
the larger fractional systematic error for Raa measure- 
ments compared to the i>2 measurements, the deviation 
of t | 2(A r par t) from the data is more dramatic than that 
for the -RAA(-Wpart)- Nevertheless, Fig. [2] clearly shows 
the importance for any model to simultaneously describe 
the Raa and the azimuthal anisotropy of the data. 



The fact that the high pt i>2 at RHIC exceeds expec- 
tation of pQCD jet-quenching models was first pointed 
out in Ref. [23| in 2002. This was not a serious issue back 
then since the px reach of early measurements was rather 
limited, and the v-i could be strongly influenced, up to 6 
GeV/c for pions, by collective flow and recombination ef- 
fects rather than jet quenching 27). Fig. [2] clearly shows 
that the v% at pr above 6 or even 9 GeV/c still exceeds 
the pQCD-based energy loss models. It is possible that 
geometrical effects due to fluctuations and CGC effects, 
ignored in these models, can increase the calculated vi\ it 
is also possible that the energy loss process in the sQGP 
has a steeper I dependence (e.g. AdS/CFT) than what is 
currently implemented in these models. 

To test whether these two ideas could bridge the differ- 
ence between data and theory, we compare the data with 
the JR model from [24[ . This model is based on a naive 
jet absorption picture with an exponential survival prob- 
ability e _K/ for jets, where the line integral I = J dl p 
is chosen for a quadratic dependence of absorption in a 
longitudinally expanding medium, and n is tuned to re- 
produce the central Raa data. The medium density p 
is given by two leading candidates of the initial geome- 
try: MC Glauber geometry p GL (x,y) = 0.43p part (x, y) + 
0.14p co n(a;, y), i.e. a mixture of participant density profile 
and binary collision profile from PHOBOS [25| : and MC 
CGC geometry pcgc( x ,v) of Dresher & Nara The 
effect of fluctuations for both profiles were included via 
the standard rotation procedure [13J. The short-dashed 



curves in Fig. [3^a) show that the result for Glauber ge- 
ometry without rotation (pen,) compares reasonably well 
with those from WHDG [22| and a version of ASW model 
from (2^|. Consequently, we use the JR model to esti- 
mate the shape distortions due to fluctuations and CGC 
effects. The results for Glauber geometry with rotation 
(PglO an d CGC geometry with rotation (pcgc) eacn ^ ea( i 
to an ~ 15 — 20% increase of V2 in mid-central collisions. 
However, these calculated results still fall below the data. 

Figure \3j[b) compares the same data with three JR 
models for the same matter profiles, but calculated for 
a line integral motivated by AdS/CFT correspondence 
I = J dl I p. The stronger I dependence for pql signifi- 
cantly increases (by > 50%) the calculated v 2 , and brings 
it close to the data for mid-central collisions. However, a 
sizable fractional difference in central bin seem to require 
additional increase from fluctuations and CGC geometry. 
Fig. [3] (b) also shows a CT model from [2(| , which im- 
plements the AdS/CFT I dependence within the ASW 
framework [29] ; it compares reasonably well with the JR 
model for pql (short-dashed curves). Note that the CT 
or JR models in Fig. [3] have been tuned independently 
to reproduce the 0-5% ir° Raa data, and they all de- 
scribe the centrality dependence of Raa very well (see 
Fig. [3] (c)-(d)). On the other hand, these models pre- 
dict a stronger suppression for dihadrons than for single 
hadrons, opposite to experimental findings 28|. thus a 
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FIG. 3: (Color online) V2 vs -/V par t in 6-9 GeV/c com- 
pared with various models: (a) WHDG 22] (shaded bands), 
ASW [H| (solid triangle), and three JR calculations [3] 
with quadratic / dependence with longitudinal expansion for 
Glauber geometry (dashed lines), rotated Glauber geom- 
etry(long dashed lines) and rotated CGC geometry (solid 
lines); (b) Same as (a) except that AdS/CFT modified calcu- 
lation in ASW framework (triangle) from [26i ] is shown and 
the JR calculations were done for cubic I dependence with 
longitudinal expansion; (c)-(d) the comparison of calculated 
i?AAS from these models with data. 

global confrontation of any model with all experimental 
observables is warranted. 

In summary, we presented results on ir° azimuthal 
anisotropy (^2) in 1 < px < 18 GeV/c in Au+Au colli- 
sions at ^/snn=200 GeV. The measurements indicate siz- 
able V2(pr) that decreases gradually for 3 < pr < 7 — 10 
GeV/c, but remains positive for px > 10 GeV/c. This 
large i>2 is striking in that it exceeds expectations of 
pQCD energy-loss models even at pr ~ 10 GeV/c. Esti- 
mates of the i>2 increase due to modifications of initial ge- 
ometry from gluon saturation effects and fluctuations in- 
dicate that they are insufficient to reconcile data and the- 
ory. Incorporating an AdS / CFT-like path- length depen- 
dence for jet quenching in a pQCD-based framework [26[ 
and a schematic model [24| both compare well with the 
data. However, more detailed study beyond these sim- 
plified models are required to quantify the nature of the 
path-length dependence. Our precision data provide key 
experimental constraints on the role of initial geometry 
and for elucidating the jet quenching mechanism. 
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